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Based on the 4D momentum, the influence of quintessence on the gravitational frequency shift and 
the deflection of light are examined in modified Schwarzschild space. We find that the frequency 
of photon depends on the state parameter of quintessence Wq: the frequency increases for — 1 < 
Wq < —1/3 and decreases for — 1/3 < < 0. Meanwhile, we adopt an integral power number a 
{a = Suq + 2) to solve the orbital equation of photon. The photon's potentials become higher with 
the decrease of uiq. The behavior of bending light depends on the state parameter uiq sensitively. In 
■ particular, for the case of Uq — —1, there is no influence on the deflection of light by quintessence. 

CNj , Else, according to the H-masers of GP-A redshift experiment and the long-baseline interferometry, 

^ ' the constraints on the quintessence field in Solar system are presented here. 

^ . 
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O^. I. INTRODUCTION 

U . 

L^j| It is well known that the universe accelerating is proved by the type la Supernovae (SNe la) [l| and the 

microwave background (CMB) Q in Wilkinson Microwave Anisotropy Probe (WMAP). The large scale distribution 
^ of galaxies 3] shows the existence of cold dark matter (CDM) with the ratio of ^Icdm = 0.27 ± 0.04. The current 
CNj universe is dominated by the dark energy component which contributes fide ~ 0.67±0.06 to the critical energy density. 
' Else, dark energy should also have a negative pressure to ensure the acceleration of universe. 

' These astronomical observations constrain current state parameter oj to be close to the cosmological constant case 
, which corresponds to a fluid with w = — 1. Actually, the observations also indicate a little bit time evolution of 

w. Thus, many people considered the corresponding situation in which the state equation of dark energy changes 
with time. From particle physics, the scalar fields is natural candidates of dark energy. A wide variety of scalar field 
dark energy models have been proposed such as quintessence [5| , K-essence tachyon field Phantom (ghost) 
. ^ [ dilatonic dark energy Q and so on. The quintessence model refers to a minimally coupled scalar field with a potential 
which decreases as the field increases. The field is coupled by the gravity through a Lagrangian taking the form 



OO 

o 



where 4> is the scalar field and V{(j)) is the potential. 

On the other hand, many people have discussed minutely how to combine the quintessence energy momentum 
tensor with Einstein equations 11 1 12 1 [3]. However, the initial solutions obtained by 11 1 12 1 have no horizon and 



no 'hair'. Hence, the black hole can not be formed in those metrics. Later, Kiselev [13| adopted a nonzero off-diagonal 
energy momentum tensor being proportional to diagonal terms, i.e. C(r) 7^ oc B{r), where C{r) and B{r) are the 
coefficients of the energy momentum tensors. The energy momentum tensor of the quintessence are 

Tl^A{r), T^^O, Ti^C{r)rl+B{r)5l. (2) 
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The constant coefficient C{r)/B{r) satisfies the additivity and hnearity conditions. So a exact black hole solution 
surrounding by quintessence matter is obtained by the static coordinates 13[ . Since then, many people have performed 
its quasinormal modes 14 1, entropy 15 1, geodesic precessionjl^ and so on. However, as far as we know, there is no 
work discussing its influence on gravity test in the view of 4D momentum. In this paper, we start from the 4D 
momentum and calculate carefully the gravitational shift and the deflection of light in the modified Schwarzschild 
space encoded the quintessence matter. From the consideration of the supernovae dimming, we should make the 
dark energy satisfies oj < —2/3 to fit the observations (IQJ . But in this paper we evaluate the quintessence matter 
with the state parameter in the range of ujq G [—1,0]. Meanwhile, we also compare the result of quintessence with 
that of cosmological constant [l3] to verify its rationality based on the fact of that ACDM is reduced from a special 
quintessence field of = — 1 to a certain extent. It should also be noted that the problems of gravitational frequency 
shift and light deflection are studied in the context of central gravitational field rather than the cosmology here. 

This paper is organized as follows: in section II, we present the Kiselev solution, i.e. the Schwarzschild black hole 
surrounded by quintessence matter. In section III, we calculate the frequency gravitational shift of photons in which 
the frequency is inverse proportional to local ^/goo■ In section IV, we study the deflection of light according to the 
different state parameter LOq. In section V, we use the astronomical observations to constrain the quintessence and 
examine whether the observable effect is big enough to be measured. Section VI is a conclusion. We adopt the 
signature (+,—,—,—) and put h, c, and G equal to unity. 



II. SCHWARZSCHILD SPACE SURROUNDED BY STATIC SPHERICALLY SYMMETRIC 

QUINTESSENCE MATTER 

Before performing the gravitational test of modified general relativity (GR), we would like to introduce Kiselev's 
black hole solution 13]. The derivation starts from a spherically symmetric static gravitational field, 

ds^ = e-^'-'-^dt^ ~ e^^'-^dr^ - r^{dd^ + sin^ 6'V)- (3) 

It assumes that the quintessence matter distributes evenly outside black hole and the energy momentum tensor takes 
the forms 

Tl = (4) 



-(l + 3i?)^+i?^^ 



(5) 



where Pq is the density of quintessence matter. The parameter B depends on the internal structure of quintessence. 
The isotropic average over the angle components is 

<Tl>=-pq(r)l^^-V^Sl, (6) 

where the relationship < r^r-' ^(5^r„r" is used and the state equation is Pq — i^qPq where uoq — 7/3. 
According to the additivity and linearity principle, i.e. = , the free parameter B{ijjq) can be fixed as 

5 = -^. (7) 

QUq 

Meanwhile, this assumption directly leads to a relationship \ -\- v = Q. Hence the energy momentum tensor can be 
written as follows, 

Tl = T; = Pq, (8) 
T'e = T^^-]^Pq{?,Uq + l). (9) 

If we set A = — ln(l + /), a differential equation of / can be obtained 

'^'0 + 3(l+'^',)'^^ + (3^9 + l)/ = O, (10) 
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which has two exact sohitions 



- (11) 



fBH = (12) 

r 

where a and Vg are the normahzation factors. Jbh is the ordinary point-hke black hole solution such as the 
Schwarzschild black hole. The quintessence density is 

So if we require the density of energy positive, pq > 0, we deduce that a is negative for ujq negative. The curvature 
has the form of 

i? = 2T,- = 3a^,i^. (14) 

Apparently, r = is the singularity for Uq ^ {—1, 0, 1/3}. Combining linearly solutions pTjl and p2p . we can get a 
general solution. 



i-^-E(v) ' (15) 



where n indicates several fields and rq is a normalization factor which has the dimension of length. The free 
quintessence creates an outer horizon of de Sitter universe at r = rg for 



l<^g<-o, (16) 



and also generates an inner horizon of black hole at r = for 

-\<OJq<Q. (17) 

Here we only consider the influence of quintessence matter on Schwarzschild space. The metric of Schwarzschild space 
is modified by a new form, 

'^^^ = (l - ^ - ^) - (l - ^ - ^) - r^de- + ^d^')^ (18) 

where M is the black hole mass. When a = this metric reduces to a pure Schwarzschild case. Otherwise, when 
LOq = —1 this metric reduces to a Schwarzschild-de Sitter one. It should be noticed that comparing with the large 
total mass of dark energy, the relatively small mass of black hole can be treated as a constant. 



III. FREQUENCY GRAVITATIONAL SHIFT 

When photon propagates in a stable gravitational field, the stationary observers at different positions obtain different 
frequencies. This phenomenology is the so called gravitational redshift which is one of the special effects purely due 
to curved space. Certainly, it is also a natural result from the equivalence principle in general relativity. For the static 
metric g^i/ (|18p . the photon's energy can be written as 

E^P^U^, (19) 

where is the 4D momentum of photon, is the 4D velocity of observer which is in the form 

= ^^(1. 0. 0: 0). (20) 
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Submitting Eq. ((20| into Eq. ((T9)) . the energy can be rewritten as 

E^^. (21) 

Using the Planck relation E =^ hv, the frequency and metric component goo yield 

^/gw)V = Pa/h. (22) 

The right hand side of Eq. (j22p is a constant since po is a conserved parameter when photon moves in this stable 
space. 

When Schwarzschild space is surrounded by the free static spherically symmetric quintessence matter, the time 
component of metric (|18p is 

2M a 

500 = 1- — -;3i:;+T- (23) 

Meanwhile, we assume that the emitter is at Xq = (ro, 0, 0, 0) and the receiver is at x = (r, 0, 0, 0). The photon 
signal has fixed spatial coordinates and their 4- velocity is tangent to the static Killing field f . Furthermore, the 
photon frequency is inverse proportion to the local Vgoo and the frequencies of emitter and receiver must satisfy 



_^ _ / .9oo(?'o) 

Submitting the lapse function Eq. (|23p into above ratio v/v{)^ we can obtain 



(24) 



a 



(25) 



^0 5^?(^o)r^'+VJ' 
where g'^Q = 1 — 2M/r and we assume 2M/r ^ a/r'^'^«+^. In the weak field limit 

M/r < 1, (26) 



Eq. can be simplified to 



where 



V . M M ^ 

— = 1 + + 27 

1^0 r ro 



- 2 L3"<J+i r-o'^'+O ■ ^^^^ 



Hence, the relative gravitational redshift (or redshift parameter) z is 

v-vq M M 



vq r ro 



(29) 



where the first and second terms are the result of GR and the last one describes the amendment to Schwarzschild 
space by quintessence matter. Furthermore, the modified term satisfies the following relation 

>0 for -l<uj„<-l/3, 

' ' (30) 
<0 for -l/3<LUg <0. 

Comparing with the usual pure Schwarzschild result we can find the redishift increases for — 1 < cj^ < — 1/3 and 
decreases for — 1/3 < < 0. 
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IV. DEFLECTION OF LIGHT 



Because that the movement of photon is described by the null line element ds = 0, s can not be treated as the affine 
parameter any more. The original definition of 4D momentum — is not suitable to describe massless photons. 

Considering the above reasons, we choose an arbitrary scalar parameter ^ as a new affine parameter. Therefore the 
4D momentum of photon is redefined as 

dC 

When photon goes in the static spherically symmetric space ([T5|) . po and are still conserved parameters. The first 
and second motion equations about t and ip are obtained naturally in the following equations: 

The third motion equation can be given by the normalization relation of photons g^uP^p'^ = 0, 

^\\e^~ — (i- — - — 



^) (34) 



The dynamic evolution of photon is determined fully by Eqs. ((5^ . ([55)1 and in this space. Combining Eq. 
with p4l) . the photon's trajectory equation can be gotten as 



1 drV f eV I [ 2M 



r'^ dip J \L J r'^\ r j-S,i^q+ij' (35) 

Meanwhile, we define two new parameters, one is the impact parameter b = L/E, which means the effective sighting 
range, the other is the photon's effective potential \/B'^{r) where 

2M - 



Bir)=r\l- — -^^,j . (36) 

The effective potential, which is illustrated by FigJU becomes higher with smaller state parameter uoq. 
According to Eq. (|36p . photons' orbital equation ([55)1 can be rewritten as 



f drV _ 1 



dp) 52 52 (r)' 



(37) 



Calculating the first order derivation of the trajectory equation (j35p with respect to p^ the photon Binet equation 
can be obtained as follows, 

j^=^u + 3.^+ '"t''t?^ .", (38) 
dip^ 2M'3"<!+i ^ ' 

where u = M/r and a = 80;^ + 2. In the right hand side of the above equation, the second term is the general relativity 
correction and the third one is the quintessence contribution. Because u can be treated as a small quantity in the 
weak field regime, we can solve Eq. (I38p by the successive approximation method. 

The negative pressure of quintessence matter {j)q — tOqPq) supplies the state parameter — 1 < cijq < 0. The 
exponential a in the additional term of Eq. ((38|) is in the scope — 1 < a < 2. It is more difficult to solve Eq. ([38|) 
just by the undetermined state parameter LOq and normalization a. In order to solve this problem, we select the 
solvable integers (a — —1, 0, 1, 2) in the range of u)q S [—1,0]. The four cases' results, which are (a — 2,ujq = 0), 
(a = 1,^5 = —1/3), (a — Q.oJq — —2/3) and (a = — l,^^ — —1) respectively, are listed in the following TABLE [J 
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FIG. 1: The effective potentials 1/B^{r) of pliotons witli a — —0.05 and unit mass M = 1. The corresponding Schwarzschild 
case a = is drawn by solid line. 
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FIG. 2: Deflection angle 5/3 versus quintessence's normalization parameter a with AI = 1, R = 10. The solid line denotes the 
pure Schwarzschild case [f 8| whose deflection angle Sp = 4GM/R corresponds to the case of a = — I. 



Here, the parameters uq, y and ?/o are uq = y = u + ^/i^SaA/^ = uq + v^IEMII^^ where R is the solar 

radius. The calculation detail are shown in Appendix [XI iBl [Cl and iDl 

The forms of influence on light deflection by quintessence are different according to various state parameters ujq. The 
deflection angles depend sensitively on quintessence's normalization parameter a. The deflection angles are illustrated 
in Fig. [2] For the case = —2/3, the deflection angle increases monotonically as \a\ increases. On the contrary, the 
deflection angles decrease with increasing \a\ in the cases of ujq ~ —1/3 and LUq = —2/3. Meanwhile, comparing to the 
pure Schwarzschild case, the deflection becomes larger in the case of a = and becomes smaller in the cases a = 1 
and 2. These variations are caused by the fact that the modified term ^'^^Sj+'i u"": in the photon's orbital equation 
psp . gives the different solutions. The behaviors of deflection heavily depend on the quantities of state parameter ujq. 
It should be notice that in order to easily analyse the results, we recover the international system of units from the 
natural unit. 
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TABLE L The influence of quintessence on deflection of light 



state parameter ujq reduced Binet equation the first order approximate solution deflection angle 

a = 2.Ljg=0 ^ + u = (3 + 1^) «^ u = uo cos¥> + (1 + 2^5-) Mo(l + sinV) = ^1^ = ^ {i + ^) 

2 ^2 

a=l,w, =-1/3 ^ + (1 - Q)ti = 3m^ u = uo cos(i^Vl - ct) + (l + sin^(¥>Vl - ")) = j^—^^ ^ 



a = O.Ljg = -2/3 j-^+u-^=3u^ i; = j/p cos[y(l - 6aM)^^''] + {1 + sin^ [y (1 - 6aAf )^/'*] } (5i^ 



— 1 , CiJq = — 1 ^ 2 -\- u — 3u^ u — uq cos + Uq ( 1 + sln^ !p) S^p — 4 



V. EXPERIMENTAL CONSTRAINTS 

We know that astronomical observations, such as the gravitational Icnsing measurements, become important in 
determining cosmological parameters. Here we will do this work in the view of black hole. We compare the theo- 
retical predictions of gravitational frequency shift and deflection of light with the typical experiments and give the 
experimental constraints on the cosmological parameters. 

Firstly, we consider the case of gravitational frequency shift. The results of frequency comparison/clock comparison 
are given by the hydrogen masers frequency standard in the GP-A redshift experiment [19i] . The continuous microwave 
signals are generated from H-masers located in the spacecraft and at an Earth station where the spacecraft was 
launched nearly vertically upward to 10000 km. This experiment reached a lO"^** accuracy. So, if we consider the 
quintessence field's effect, the parameter a has to satisfy the constraint 

|a|<2x 10-14 X ^ °^ , . (39) 

The constraints of field parameter a are presented in TABLE [ill Generally speaking, the upper limit of field 
parameter \a\ increases with bigger state parameter LOq. When ujq crosses over —1 the maximum |a| is below the value 
of 10~^^, which is consistent with the cosmological constant case \a\ < 10~^^ [l3]- However, there is a singular point 
located at = ^1/3 since the definition of ujq — —1/3 does not exist in Eq. ipS)) . Else, when Wq —1/3, i.e. A^, 0, 
the redshift parameter z ([29|l reduces to the result of pure Schwarzschild case. 

TABLE II: The constraint on field parameter a by H-masers of GP-A redshift experiment. 

ujg -0.2 -0.4 -2/3 -0.8 -1 

Estimate on Q |q|<3x10"'' |a|<6xl0"" |q|<8x10~^'^ |q| < 5 x 10"^^ |q| < 6 x 10"^'' |q| < 3 x 10^^** 



Secondly, we consider the other case of deflection of light. In order to obtain experimental constraints from the 
light deflection result, the final light deflection angle Sip should be expressed in terms of the deviation A^jj from the 
general relativity prediction 5lpgr for the Sun. 

5ip = 5^gr{^ + Ald), (40) 

where SipcR — AGM/R. The best available constraints on SipcR. come from long-baseline radio interferometry which 
shows that |Ald| < 0.0017 "20]. Submitting the deflection angles of ujq = -1, -2/3, -1/3, shown in TABLE |I] 
into Eq. (|40|) . we can obtain the corresponding constraints on the field parameter a, which are present in TABLE 
mil The constraint on a becomes stronger with the decreasing state parameter iVq. 



VI. CONCLUSION 



In this paper, we have studied the influence of free quin-tessence on gravitational frequency shift and deflection of 
light based on the 4D momentum. We summarize what have been achieved. 
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TABLE III: Estimates on a from Solar system observation. 



Uq A_Li5 Estimate on a 

-1/3 1-(1-q)-3/2 |a| < 10-3 

-2/3 (1 + f aM)(l - ^) - 1 |a| < 10"^^ 
-1 — 



1. The influence of quintessence matter on our space is expressed mathematically as a new metric form. Taking 
advantage of the symmetries of Kiselev solution, we can avoided to solve directly the original geodesic equation. The 
inner product M°^a is constant along the geodesic in a Killing field with a geodesic tangent u"" . In this paper, 
starting from the 4D momentum the geodesies problem is reduced to the problem of one-dimensional motion of a 
particle in an effective potential. The n ull g eodesic of the modified Schwarzschild geometry containing the quintessence 
matter is solved by the usual method We analyze the behavior of light ray in the weak field regime r ^ M. 
Also, this result can be applied to the exterior field of ordinary body such as the Sun or more large astronomical scale. 
Meanwhile, the quintessence matter is also treated as a very small component in whole space, i.e. M/r ^ a/r^'^'''^^ . 
These appropriate approximation can help us to simplify the calculation of the gravitational redshift and the bending 
of fight. 

2. The small goo in a metric is corresponding to the high frequency because the frequency of photon is inverse 
proportional to local y/goa, i.e. cx l/^.9oo- The gravitational redshift expression ()29|1 is obtained naturally in the 
weak field Hmit. It is clearly that the frequency is larger (or smaller) for —1 < ojq < —1/3 (or — 1/3 < < 0) than 
that of the pure Schwarzschild case from the result of redshif (|27|) . 

3. Comparing with the gravitational redshift, the deflection of light is more or less complex with the uncertain power 
a in the orbital equation ([55]) . We choose the solvable integral numbers a in the range of e [— 1 , 0] to obtain the 
deflected angle. The solutions of the orbital equation ([55]) depend sensitively on the value of usg. Furthermore, when 
a = — 1, i.e. LOq — —1, there is no influence on bending of light by quintessence matter, which is in accordance with 
the theoretical prediction of SdS case [l3| • This point also illustrates well that the border case of the extraordinary 
quintessence ujq = —1 covers the cosmological constant term 

4. The parameter a is of course a crucial parameter, which indicates the influence of quintessence matter over the 
space, it is therefore necessary to give its realistic order of magnitude. Based on the the relation (fT^]) . we evaluate 
roughly the order of a in TABLE llVl in the four different astronomical scales, i.e. Solar, Galaxy, Cluster of Galaxies 
and Supercluster. In the range of Solar system, the measuring unit of the distance between two planets is Astronomical 
Unit (AU). It is well known that the nearest Mercury is about 0.4 AU and the most remote Pluto is about 40 AU far 
from the sun. Here the average orbital radius of Pluto is chosen as the value of the parameter r in Solar system. The 
other three large astronomical scales are estimated in the order of magnitude. If quintessence matter can be treated 
as the dark energy model, its density pq can be assumed as the critical density of the universe, namely, 

{10~46 GeV^ {natural unit), 
10^"* eV/cm^ {energy density unit), (41) 
10^^^ g/am? {mass density unit). 

Hence, substituting the density pq and the different scales r into relationship (|13[) . we can get the order of magnitude 
of a with different state parameter LOq, which is shown in TABLE IIVI Apparently, the quintessence matter makes 
more influence in the large scale. Comparing the long-baseline radio interferometry experiment results TABLE IIIII 
with the theoretical expectation values TABLE IIV[ we can find that the case of ujq = —2/3 violates many orders of 
magnitude of \a\ in solar. So the case of LOq — —2/3 should be abandoned in solar system. 



9 



5 There are two easily confused points should be clarified in the end. Firstly, the parameter a in Ea. ((TT|) and 
Eq. p^ are the same as a matter of fact. The original Kiselev's black hole solution is the general form of exact 
spherically-symmetric solutions for the Einstein equations describing black holes surrounded by the quintessential 
matter with the energy momentum tensor yielding the additive and linear conditions ([8]) and ([9]), i.e. Eqs.(13) — (14) 



in original paper 13|. If we only consider the simplest singular Schwarzschild case, the general solution to Eq.([T 



should be in the form of 

e-^ = l-^--^- (42) 



r 



So the modified Schwarzschild black hole solution (fT8|) can be obtained without changing parameter a. Mathematically, 
if we adopt a singular field, not the mulriple fields, the first term of Kiselev's solution should be 



Comparing with the exact solutions (|lip and p2)) . we can get following equation 

So it is clear that two relationships r^"^"^^ — a and ujq — ujq exist here. Hence, the parameters a in Egs. pTj) and ([18]) 
are the same. 

Secondly, the parameter M is actually the mass of black hole (or the center mass of gravitational field) , not the 
mass of outside quintessence matter. The form of quin-tessence field is denoted by the energy momentum tensor 
yielding the conditions (8) and (9). Since the dark energy almost has the ratio Q.de — 0.67 ±0.06 in universe, the mass 
of quintessence Mq is far more than the mass of a stellar black hole M, i.e. Mq ^ M . So the quintessence density p 
(fT3|) does not depend largely on the mass of black hole and M can be assumed as a constant. This assumption also 



can be found in the corresponding extensions of this quintessence black hole solution 1J|, 15|, [1 



TABLE IV: The theoretical magnitude order of parameter a in several astronomical scales 
astronomical scales Solar (40AU) Milky Way (lO^Pc) Cluster of Galaxies (IMPc) Supercluster (lOMPc) 



^5 



= -1/3 7.16 X 10"^^ Ag.m~M.90 X 10^ Kg.m-^ 1.90 x 10^^ Ap . 1.90 x 10^'' Ap . 

= -2/3 5.98 X 10-23 Kg . m,-^ 3.08 x lO'^^Kg . -m.-^ 3.08 x 10"^" Ag . m'^ 3.08 x 10~" A3 . m' 

= -1 6.67 X 10-33 Kg.m-^ 
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When the state parameter = 0, Eq. pS]) reduces to 

Because of the small m, we first ignore the 2th order small quantity iv? . The zeroth order approximate solution is 

u = uq cos if, (A-2) 
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where uq = GM/R and R is the solar radius. This is a straight hne normal to pole axis. Substituting Eq. (|A2|1 into 
the right hand side of Eq. (|A1[) . we can get 

:73 + " = 3 + ^ ho cos' ip. (A3) 



This equation has a particular solution, 



l + ^)"o(l + sinV)- (A4) 



So the first order approximate solution to Eq. (jAip is 

u = itocos.^ + (^1 + Uo(l + sin^). (A5) 

The azimuth angle of null order approximate solutions (|A2p are ±7r/2 in very far place u = Q. However, the azimuth 
angles of the first order one (|A5p are ±(7r/2 + /3) in u = 0. The deflection angle /? is a small quantity which satisfies 



Mo 



sin/3+ (1+2^)^^(1 + cos2/3) = 0. (A6) 



We use the Taylor expansions of sine and cosine functions and keep the basic term. The deflection angle is written 
formally as 

Comparing with the result of pure Schwarzschild space [l^, the additional term containing parameter a is the 
contribution of quintessence matter. The case of = affects the light deflection only through the particular 
solution (|A4p . But the zeroth order approximate (|A2|) is unchanged. 

APPENDIX B: THE CASE OF a = l,^;, = -1/3 
When = -1/3, Eq. ([38]) reduces to 

— + 1 - a u = ?,v?. Bl 
The process of solving this equation is similar to the formal one. The zeroth order approximate solution is 



u = Uq coa{(p\/l — a). (B2) 

The particular solution is 



u=-^(l + sm''{ipVT^)). (B3) 
1 — a 

The first order approximate solution is 



u = uq cos(^V1 - a) H — (1 + sin' {ipVl - a)) . (B4) 

1 — a 

So the deflection angle is 

A 4 GM 

The deflection behavior of the case a;^ = — 1/3 is different from the case of = 0. On the contrary, the quintessence 
matter effects the photons deflection through the zero order solution (|B2p . So the later particular solution and first 
order approximate solution vary accordingly. 
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APPENDIX C: THE CASE OF a = 0,a;, = -2/3 
When = -2/3, Eq. ([38]) reduces to 

^+u- — = ,u^. (CI) 

In this case, the additional term aM /2 changes the form of the zeroth order approximate solution. In order to simplify 
calculation, we use a transformation 

1 - VI - 6aM 

u = y+ . (C2) 

So the Eq. (jCip can be rewritten as a usual orbital equation 

Its first order approximate solution is 



Vl - 6aMy = Sy^. (C3) 



y = 2/oCOs[(p(l - 6aM)i/4] + ^° |l + sin2[^(l - 6aM)i/4]\ (C4) 

Here we only consider the small contribution of quintessence matter in the large space background. Hence, when 
u — > 0, the limit of y is also small. So the deflection angle is 

'^^ ^ (l-6aM)3/^ ' 

where 



VI - 6aAf - 1 

yo = uo H . (C6) 



The additional term vanishes for the factor LUq + 1 = 0. Hence, the influence of quintessence can be ignored. This 
point is also justified by potential analysis. Since the exponential of additional term is negative Sujq + 2 < 0, should 
not be considered as a small quantity any more. The successive approximation is failing here. However, the photon's 
effective potential can be rewritten as 



1 



-2 I 2M 2 
= r II r a 



The photons' trajectory equation ([35]) of ujq = —1 is 



dr 

dip 



52 B^{r) 



1/2 



(Dl) 



(D2) 



where b is the impact parameter and the ± denotes to increasing or decreasing r. According to the condition of 
perihelion position tq, we can get 



dr 
dip 



1 



1 



0. 



(D3) 



So the impact parameter b can be expressed by perihelion position tq. Submitting b into Eq. (|D2p . the quintessence 
matter can be dropped out, 

1/2 

(D4) 



dr ^2 

dip 



111 

ro 



2M 



1 2M 
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Therefore, the quintessence matter has no influence on the hght deflection in this case, which is also justified by the 
articles [17 1. 
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